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Abstract 

We use particle image velocimetry to measure the sedimentation dynamics of a semi-dilute 
suspension of non-Brownian spheres at Reynolds numbers, 0.001 < Re < 2.3, extending from the 
Stokes to the moderately inertial regime. We find that the onset of inertial corrections to Stokes 
sedimentation occurs when the inertial screening length I = a /Re becomes similar to the Stokes 
sedimentation length ^o, at Rcc = a/^o 0.05. For Re > Rcc, inertial screening significantly 
reduces both the magnitude and spatial extent of the particle velocity fiuctuations. A modified 
Hinch force balance model connects the fluctuation magnitudes cry/V to the correlation sizes 4. 



The sedimentation dynamics of non-Brownian spheres is a fundamental problem in 
physics [1] and has been the subject of intense activity in recent years (for a review see [2]). 
Much of the focus has been on particles slowly settling in viscous liquids, conditions that 
correspond to very low Reynolds numbers, Re <^ 1, called the Stokes regime. The Reynolds 
number is the ratio of inertial to viscous forces in fluids, and is defined as Re = 2Va/ v {V is 
the particle velocity, a the radius, and v the fluid kinematic viscosity). In the Stokes regime, 
where inertial forces are insignificant, concentrations of spheres display i?e-independent large 
amplitude velocity fluctuations oy jV during settling [Sj HI El El [7]. Significantly, experi- 
ments O IHl [71 [8] , simulations [9], and theory [10] have shown that the fluctuations display a 
characteristic spatial size ^, despite the fact that the hydrodynamic interactions emanating 
from a single isolated sphere are of infinite range (oc 1/r) [11]. The Stokes screening length 
^ represents regions of concentration fluctuations, a^, that drive velocity fluctuations cxy 
and particle diffusion D ~ cry^. While an explanation of the origin of screening in Stokes 
sedimentation remains controversial, its existence and central importance for a description 
of sedimentation is not. 

The sedimentation dynamics at higher speed flows where inertial forces become signif- 
icant. Re ~ 1, have received much less attention despite its fundamental importance and 
widespread relevance to numerous chemical industries ^[2] . In contrast to the Stokes regime, 
for an isolated sphere falling at moderate Re, there is a distance beyond which the 1/r Stokes 
like hydrodynamics are screened [TT]. The inertial screening length, / ~ a/ Re, becomes 
shorter as the Reynolds number is increased. For concentrations of settling spheres, the 
Reynolds number dependence of the hydrodynamic interactions brought about by inertia is 
predicted to result in Reynolds number dependent sedimentation dynamics as well [13]. Ex- 
periments designed to investigate these changes, however, are greatly lacking. We are aware 
of only a single experiment, by Cowan et al. [H] using novel ultrasonic techniques, that 
has examined the fluctuation dynamics of spheres at moderate Re ^ 1. Surprisingly, they 
concluded that the velocity fluctuations were independent of Reynolds number for Re ^ 1. 
In the absence of experiments describing the sedimentation dynamics beyond the Stokes 
regime, our understanding of the basic physics, and our ability to develop and test model 
theories, remains very limited. 

In this letter, we describe experiments that demonstrate how moderate amounts of fluid 
inertia can significantly affect the settling dynamics of spheres. The onset of inertial correc- 



2 



tions to Stokes sedimentation occurs when the inertial screening length / = a /Re becomes 
as small as the Stokes sedimentation length ,^0; at Rcc = a/^o- For > R^c, inertial 
screening significantly reduces both the magnitude and spatial extent of the particle veloc- 
ity fiuctuations. A modified Hinch force balance model, with Reynolds number dependent 
drag coefficients, connects the fluctuation magnitudes cry/V to the correlation sizes ^ over 
the entire range studied, 0.001 < Re < 2.3. 

The particles used in our experiments were monodisperse glass beads of radius a = 137±9 
lim. They were dispersed in various mixtures of glycerol and water, enabling us to examine 
the very low to moderate Reynolds number regimes, 0.001 < Re < 2.3. In all cases the 
volume fraction is = 0.06. The sample cell was a rectangular glass tube of dimension 8 x 
80 X 305 mm, and the temperature was at the ambient value T = 23±1'^C. Particle velocities 
were measured using a particle image velocimetry (PIV) [15J apparatus consisting of some 
specialized image processing software and hardware purchased from Dantec Instruments. A 
large cross section of the cell was imaged (3x3 cm), so that several thousand particles could 
be simultaneously studied. Initially, random dispersions were prepared by vigorous shaking 
of the cells. The location of the imaging window was far from both the sedimentation front 
and the sediment growth. Each velocity fields is a map of 41 x 42 vectors, and each velocity 
vector is the local average of two to four spheres. 

We begin by examining typical particle velocity fields during sedimentation. Figures 
[T]^a-c) show results for particle velocities Vj from three samples spanning the very low. 
Re = 0.001 <^ 1, to moderate. Re = 2.3, Reynolds number regimes. In the Stokes fiow 
regime. Fig. [T]^a), the patterns look similar to those reported in the literature O [7], with 
large magnitude fiuctuations occurring in extended swirling regions. In the presence of some 
inertia. Re = 0.2 in Fig. [ij^b), the pattern appears slightly more uniform, indicating a slight 
reduction in the fiuctuations relative to that seen at Re = 0.001. At the highest fiow speed. 
Re = 2.3 in Fig. ^c), there is a dramatic change in the dynamics. The fiuctuations are 
largely absent, the particles appear to be nearly falling with the same speed. We interpret 
these changes as evidence of an inertial screening of the velocity fiuctuations. 

To examine the fiuctuations in more detail. Figs. [l]^d-f) shows the same set of maps 
with the mean velocities subtracted off, 6\i = Vj — (Vj). All of the patterns in fact show 
qualitatively similar swirling regions. The sizes of the regions appear to decrease with Re. 
The dramatic decrease in fluctuations seen in Fig. flVc) is reflected in the greatly magnifled 
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FIG. 1: PIV results for three suspensions of sedimenting spheres at cj) = 0.06 and Reynolds numbers 
(a,d) Re = 0.001, (b,e) Re = 0.2, and (c,f) Re = 2.3. (a-c) Velocities Vj. (d-f) Velocity fluctuations 
calculated from (a-c) using (5Vi = Vj — (Vi). (Note the magnified scales in (e,f)). The single vector 
to the right of each map represents the mean velocity (Vj). 

scale (7 x SVi) of Fig. |l[f) relative to Fig. |ljd). 

By collecting large numbers of velocity maps over the majority of each falling column, 
we can access ensemble averaged information. We first examine the spatial correlations of 
the velocity fluctuations. The normalized autocorrelation function of the z component (|| to 
gravity) of the velocity fluctuations is defined as Cz{r) = {SVz{0)6Vz{r))/ {6Vz{0)'^), where 
(...) represents an ensemble average of several hundred vector maps. The distance vector r 
is either in the direction parallel to gravity, Cz{z), or perpendicular to it, Cz{x). 

To search for evidence of inertial screening in the spatial velocity correlations. Fig. [2] 
compares the correlation functions between a sample in the Stokes regime. Re = 0.001, 
and one with significant inertia. Re = 2.3. In all cases the functions decay to zero at large 
distances, indicative of finite range correlations. The decay lengths, however, are much 
reduced in the inertial sample. The form of the perpendicular correlation function is also 
affected, the long-range negative correlation dip in Cz{x) is completely absent in the inertial 
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FIG. 2: Spatial correlation functions of the z component of the velocity fluctuations as a function 
of distance ||, Cz{z), and _L, C^(x), to the z— axis for Re = 0.001 and Re = 2.3. 

case. 

To quantify the range of the velocity correlation functions, we fit to the empirical forms 
Cz{z) = exp[—{z/^y-^] and Cz{x) = exp[— as shown in Fig. [2} Results are given 
in Fig. for ^ and from all of our samples ranging from 0.001 < Re < 2.3. The 

longitudinal length ^, the longer of the two, shows the most change with Re. For Re ^ 0.05, 
^ is independent of Re as expected in the Stokes regime, and the values are all in good 
agreement with the scaling relation ~ lla^"^/^ previously found in very low Reynolds 
number sedimentation [5]. For Re ^ 0.05, the behavior begins to change, and ^ follows a 
roughly logarithmic decay with Re. The decrease of ^ with Re demonstrates that inertial 
screening occurs at shorter distances for higher Reynolds numbers. In agreement with this, 
the short range transverse correlations don't show any inertial influence until much higher 
Reynolds numbers. Re = 2.3. For analysis purposes below, we note that ^ can be well fit by 

W(l + 5.5i?e)i/^ (1) 

The velocity maps in Fig. [T] also show dramatic changes to the magnitudes of the 
velocity fluctuations as Re is increased. To quantify this we consider the ensemble averaged 
rms velocity fluctuations, = ^ {[Vi^z — V^), with results for the normalized values a^/V 
shown in Fig. ^h). oy jV is independent of Re for Re ^ 0.05, as expected in the Stokes 
regime, and the value oy jV ~ 0.75 is in good agreement with results found in Stokes flow 
sedimentation [U El [6]. In a similar way to that seen for ^ vs. -Re, the onset of inertial 
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FIG. 3: Inertial screening of velocity fluctuations, (a) Spatial correlation lengths and as a 
function of Reynolds number Re. The line through ^ is the empirical fit Eq. ([T]). The line through 
is a guide to the eye. (b) Normalized velocity fiuctuations cTy/y vs. Re. The solid line is the 
Hinch model Eq. ([2]). The dashed line is Eq. ^ with Re = 0. Inset: ^ and {ay/V), normalized 
by their Stokes values, vs. (^/a)i?e. 

influence occurs at Re ~ 0.05, and at higher speeds the fluctuations sharply decline with 
Re. At our highest speed. Re = 2.3, <Jy/V is reduced by ^ 80%. 

We now address whether the dependencies found for the correlation lengths and the 
fluctuation magnitudes can be interrelated by a simple extension of the Hinch model [U [16] 
that was originally developed for the Stokes regime [H El |6]. Within this model, velocity 
fluctuations arise from particle density fluctuations in regions whose spatial extent is the 
correlation length ^. For a random particle conflguration, the average number of particles in 
a region of size is = C,^(f)/vp, where Vp = Ana^/S. The fluctuations in number are AA^^^ = 
^A^, and mass Am^ = AN^Vp{ppartide — P fluid)- In steady state, the buoyancy force acting 
on these regions, Fg = Am^g, is equal to the viscous drag force Fd = 6n{l + j3Re)vi'^^ ^ 
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yielding AV^ = Am^g /QTiil + (3Re)fl^- The term j3Re ~ is the inertial part of the 

drag force, vahd for Re ^ 30 [I7j. Additionally, because the drag force applies to regions of 
size ^, the appropriate Reynolds number is Re^ = 2V^/v = {^/a)Re. The modified Hinch 
model then becomes 

av/V ^ 0.6[1 + O.133([e/a]i?e)°-'«]-iy^e07a, (2) 

where the bracketed term represents the influence of particle inertia, and the Hinch model is 
recovered when Re ^ [B]. (We also neglect a prefactor 7(0) = S{4>, 0) 1.0 

in our semi-dilute samples). 

To test this model, we input our fit expression for ^ from Eq. ([!]) into Eq. ^ and directly 
compare the predicted values for (Jy/V with our data. As seen in Fig. |3]^b), the agreement 
over the entire range of Reynolds numbers is remarkably good, particularly considering 
the relative simplicity of the model. The agreement shows that the drop in fluctuation 
magnitudes with Re is due to both (i) the decreasing correlation lengths ^ and (ii) the 
increasing drag force term Pr^. To illustrate their relative importance. Fig. |3]^b) shows the 
model predictions when the inertial drag force corrections are neglected {Re = 0, dashed 
line). The result shows only a slight decrease with Re, now due solely to the reduction in ^, 
and greatly underestimates the observed drop in fluctuation magnitudes. This shows that 
both terms (i) and (ii) are of similar importance in capturing the behavior of fluctuations 
at moderate Reynolds numbers. 

The value of the Reynolds number. Re ~ 0.05, at which inertial effects first appear can 
be rationalized in two different ways. First, it seems reasonable to expect that when the 
inertial screening length / becomes similar to the Stokes correlation length ^o, the fluctuations 
will show an inertial influence. This translates to an onset value I = a/Rec ~ ^0, or 
Rec ~ a/^o ~ 0.04, in excellent agreement with our estimate of Re ~ 0.05 from Fig. [3] 

Interestingly, Brenner proposed a different criterion based solely upon a comparison of 
the particle diffusion coefficient D and the solution viscosity u [18j. Particle diffusion is 
driven by velocity fluctuations, and can be estimated from D = 0.4^(jy [U El [19]. In the 
viscous Stokes regime D -C z/. Brenner argued that if the viscosity is reduced to the point 
where particle diffusion becomes similar to (viscous) momentum diffusion, D ~ i/, particles 
will diffuse faster than the momentum they are releasing into the fluid, so that the purely 
viscous Stokes flow conditions no longer apply. To test this criterion, we plot in Fig. |4]our 
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FIG. 4: Fluid kinematic viscosity particle diffusion coefficient D = OAay^, and their ratio v/D 
(the Schmidt number Sc) vs. Re. The sedimentation dynamics change when D becomes as large 
as u. The solid lines are a guide to the eye. 

values for u, D, and the ratio u/D. The results generally confirm the Brenner picture. The 
region where z//D ~ 1 is indeed where significant deviations to Stokes behavior are seen. 

Finally, to demonstrate the generality of our results, we re-examine the experiments and 
conclusions of Cowan et al. [H]. They compared many concentrations 0.19 < < 0.5, at 
both Reo = 0.007 and Reo = 0.3, and concluded that, contrary to expectations, the velocity 
fluctuations were independent of Reynolds number for Re < 1. We first note that they 
defined Rcq = 2aVo/?7o based upon the infinite dilution values Vq and rjo, not the values at 
the high concentrations studied. Using our definition Re = 2aV/ri, we estimate that the 
Rcq = 0.3 samples range from 0.001 ^ Re ^ 0.05. Our critical Reynolds number for the 
appearance of inertial influence, Rec = ci/C,o, can also be evaluated from their findings that 
^ ~ lla0~^/^, resulting in Rec ~ 0.06. We therefore find that Re ^ Rec for all of their 
samples, consistent with their findings of no inertial effects.. 

The results described here show that the onset of inertial corrections to Stokes sedimen- 
tation occurs when the inertial screening length / = a /Re becomes as small as the Stokes 
velocity correlation length ,^0) at Reynolds number Rec = a/C,o- At higher Reynolds num- 
bers, inertial screening reduces the size and magnitude of the velocity correlations. The 
success of the Hinch model in connecting cry/V with ^ suggests that the models underlying 
assumption of a random particle density distribution remains valid in the presence of inertia. 
These results provide an important benchmark for future theoretical work on moderately 
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inertial systems. 

We thank Tony Ladd for many fruitful discussions. 
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